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ABSTRACT 

We calculate the unseen-but-true GRB-event rate from a current flux-limited 
sample of 33 GRBs with individually measured redshifts. We consider a GRB- 
event rate which is proportional to the star-formation rate, in view of the GRB- 
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SNe association in SN1998bw and GRB030329. By fitting a log-normal distribu- 
^ tion of GRB peak-luminosities, we find a ratio l/f r ~ 450 for the true-to-observed 

GRB-event rates. This provides an independent derivation of the GRB-beaming 
factor 1/fb ~ 500 obtained by Frail et al. (2001) from sources with standard 
GRB-energies. We discuss applications to GRB980425. 

Subject headings: gamma-rays: bursts, beaming 
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•rH . 1. Introduction 

X 

New clues are emerging that long GRB are associated with supernovae (Galama et al. 
1998; Bloom et al. 1999; Kulkarni et al. 2000; Reichart 2001). This is illustrated most recently 
by optical emissions lines in the late-time light-curve of the HETE-II burst GRB 030329 
(Stanek et al. 2003), which are remarkably similar to those observed in GRB980425/SN1998bw 
(Galama et al. 1998). A GRB-supernova assocation provides important support for the col- 
lapsar model of GRBs, representing a violent death of evolved massive stars (Woosley 1993; 
Paccynski 1998). The short lifespan of tens of Myrs of massive stars implies that GRBs take 
place in star-forming regions (Paccynski 1998; Fruchter et al. 1999), and hence more broadly 
points towards an association to molecular clouds. The event rate of GRBs per unit cosmo- 
logical volume is hereby expected to be correlated to the cosmic star-formation rate (e.g., 
Blain & Natarajan (2000); Berger et al. (2003); Choudhury & Srianand (2002)). Because 
observations by past and current experiments (see Table I) are flux-limited, the observed 
GRB-event rate is strongly biased towards events at lower redshifts. 
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In this letter, we calculate the ratio l/f r of unseen-to-observed GRBs. Our analysis 
is based on the bias in the observed GRB-event rate towards low redshift events in a flux- 
limited sample of 33 GRBs with individually measured redshifts. We propose to use a linear 
relationship for the correlation between the true-but-unseen GRB event rate and the cosmic 
SFR, to derive from this bias a best-fit log-normal GRB peak-luminosity function. While 
it is an open question to what the degree cosmological evolution of metallicity affects the 
SFR independently of the GRB event rate (see, e.g., Heger et al. (2002)), linear relationship 
between the two will serve as a leading-order approximation. 

In §2, we discuss the cosmological model. §3 tabulates the current sample of 33 GRBs 
with individually measured redshifts. In §4, we present the best fit of the data for a log- 
normal distribution function, estimate the fraction of unseen-to-observed GRBs, and com- 
pare our results with the beaming factor 1/ f b = 500 of Frail et al.(2001), based on a sample 
of GRBs with achromatic breaks in their light-curves. We summarize our results and draw 
our main conclusions in §5. 



2. The cosmological model 

We model the intrinsic GRB-event rate in a flat A-dominated cold dark matter cosmol- 
ogy with closure energy densities £7a = 0.70 and VL m = 0.30. These values are suggested by 
BOOMERANG and MAXIMA (de Bernardis et al. 2000; Hanany et al. 2000) on the power 
spectra of the CMB and distant Type la SNe (Perlmutter et al. 1999; Schmidt et al. 1998). 
The Hubble parameter H c is taken to be 73kms _1 Mpc _1 (Freedman et al., 2001). 

Porciani & Madau (2001) provide three models of the cosmic SFR up to redshifts z ~ 5, 
reflecting some uncertainties in SFR estimates. No significant changes have been noticed 
in application of these three models, and we present only results derived from the second 
model, SFR2. We have 

R c (z) = iW*) ( f M yr" 1 Mpc" 3 (1) 

by transformation of the Rsf2 in a matter-dominated unverse (Q m = 1), 

R SF2(z) = (2 Q 2 1 ^ e e 3.t ) M yr" 1 M P C ~ 3 > (2) 

where E(Q h z) = [(1 + zf{l + zVt m ) - z(2 + z)^} 1 / 2 . 

For a GRB-event rate locked to the SFR, the true event rate between z and z + dz, as 
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observed in the observer's frame of reference, satisfies 



dRcRB{z) = Xgrb 



R c {z) dV 



dz 



(3) 



1 + z dz 



where Xgrb is the formation mass fraction of the source progenitors. Here, the division by 
1 + z accounts for time-dilatation by cosmic expansion. The element of comoving volume is 



The GRB redshift probability density can be written as (Coward et al. 2002) p(z) = 
dRcRB/dz {^j^ dRcRB / dz dz^j . Assuming that the mass fraction of GRBs progenitors is 
redshift independent, the scaling factor Xgrb is the only free parameter of our model. For 
the flux-limited experiments listed in Table I we define a probability-density function of de- 
tection as a function of redshift as Pdetect{ z ) = dRd e tect/dz ^ J Q 5 dR detect / dz dz^j , where the 
dependence on the luminosity has been integrated out in the detected event rate (Bromm & 
Loeb 2002) 



Here, p(L) refers to the intrinsic GRBs luminosity function in the BATSE energy range 
50 - 300 keV. 

The luminosity threshold as a function of redshift is given by Lu m (z) = Aird\(z)Su m , 
where d^ is the luminosity distance to a source at redshift z and where Su m denotes the 
sensitivity threshold of the instrument. Following (Bromm & Loeb 2002), we take a flux- 
density threshold of BASTE of 0.2 photon cm -2 s _1 . 



As of the time of this writing, we have 33 GRBs with individually measured redshifts 
(Table I). They collectively represent a variety of past and current experiments. Due to 
a flux threshold in each of the instruments used, the observed GRB redshift distribution is 
strongly biased towards low redshifts. This introduces a quantitative selection effect, relative 
to the redshift distribution predicted by the SFR model (Fig. 1) - the true distribution of 
GRB-redshifts representing what would be observed in the ideal case of a zero-flux threshold 
in the instruments. In view of the various instruments involved, the dependency of the 
observed GRB redhift-distribution function on the flux-limit is simulated in Fig. 2. 
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3. A redshift sample of 33 GRBs 
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4. A fit to a log-normal distribution function 

We shall assume that the GRB-luminosity function is redshift independent, i.e., without 
cosmological evolution of the nature of its progenitors. We take a log-normal probability 
density for the luminosity shape-function, with mean fx and width a given by 

where L is normalized with respect to lcm _2 s _1 . The optimized parameters of our model 
are 

(fx,a) = (124, 3) ±(2, -0.4). (7) 

This notation means that the estimated parameters can be either (122,3.4), (123,3.2), (125,2.8), 
or (126,2.6), but not (122,2.6) for instance. Our results are compatible to the expectations 
of Sethi et Bhargavi (2001) who derive a log- normal luminosity function with fx = 129 and 
o" = 2 from a different flux-limited sample. They mentioned that because of selection effects, 
the inferred average luminosity over-estimates the true mean luminosity by a factor of 2 or 3, 
and the variance is within 40-50% of the true variance. The observed and predicted redshift 
distribution, based on the SFR, are shown in Fig. 1 in case of optimal parameters (7). The 
results indicate a good fit to the data, suggesting that selection effects are adequaly modeled. 
The fraction of detectable GRBs as a function of redshift 

shows a the steep decrease as the luminosity threshold increases, making high-redshift GRBs 
less likely to be detected. 

A satisfactory test for our model is provided by comparing the predicted flux distribution 
of GRBs above the sensitivity threshold with the peak flux distribution of a sample of 67 
GRBs observed with IPN (Fig. 3). The fluxes derived from our luminosity function in 50- 
300 keV have been extrapolated to the IPN range of 25-100 keV, assuming an E~ 2 energy 
spectrum and using the formula given in Appendix B of Sethi et Bhargavi (2001). The 
conversion factor from ergcm _2 s _1 to photon cm -2 s _1 has been taken to be 0.87 x 10~ 7 , 
and the sensitivity threshold equal to 5 photon cm -2 s _1 (Hurley, private communication). 

A quantitative result provided by our model is the ratio 1/ f r of the true GRB-event rate 
(what would be seen in case of a zero-flux limit in the detector) and the observed GRB-event 
rate. For the optimal parameters (7), we have 

1/fr = 450. (9) 
The factor l/f r is between 200-1200 in the error box of (7). 
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5. A correlation between l/f r and l/f b 

Our estimate of the true-to-observed GRB-event rate l/f r is strikingly similar to the 
GRB-beaming factor l/fb of about 500 derived by Frail et al. (2001). Our analysis is in- 
dependent of the mechanism providing a broad distribution in GRB luminosities. Without 
further input, our results may reflect (a) isotropic sources with greatly varying energy out- 
put, (b) beamed sources with standard energy output and varying opening angles, or (c) 
anisotropic but geometrically standard sources (Rossi et al. 2002; Zhang & Meszaros 2002). 

We find that the GRB peak-luminosities and beaming are correlated. To see this, we 
simply note that the case of no correlation between peak-luminosities and beaming give rise 
to an unseen-but-true GRB event rate which is 1/ f r x 1/ f b ~ 2.5 x 10 5 times the observed 
rate. The true GRB-event rate hereby approaches that of Type II supernovae - we discard 
this possibility. A correlation between peak-luminosities and beaming is naturally expected 
sources (b) and (c) with standard energy output - the picture that bears out of Frail et 
al. (2001). This introduces a correlation L oc between peak-luminosity and beaming 
factor in (b) and, for a flux-limited sample, also in (c). For a flux-limited sample, both 
(b) and (c) give rise to an anticorrelation between inferred beaming and distance such that 
to leading-order 9jZ ~const. Fig. 4 shows that this anticorrelation holds in the sample of 
Frail et al. (2001). For a related discussion on estimating the GRB beaming factor from flux 
limited surveys, see Levinson et al. (2002). 



6. Conclusions 

A table of 33 GRBs with individually determined redshifts allows an estimate of the 
GRB-luminosity function, based on constant of proportionality between the GRB event rate 
and the cosmic star-formation rate. We have tested our fit by reproducing the distribution 
of peak luminosities in the IPN sample of 67 GRBs. 

A flux-limited sample introduces a ratio 1/ f r of unseen GRBs, whose emissions fall below 
the detector threshold, to observed GRBs. A best fit analysis of the luminosity function gives 
l/f r ~ 450. This number is both large and close to the GRB beaming factor of 500. The 
possibility that 1/ f r and l/fb are uncorrelated gives rise to prohibitively large true GRB 
event rates, and is discarded. The results support a relation L oc l/fb between the peak- 
luminosity and the beaming factor in the context of conical or anisotropic GRB-emissions 
with standard GRB-energies. This points towards an anticorrelation between distance and 
half-opening angle, which is approximately supported by the sample of Frail et al. (2001). 

The results presented here and those of Frail et al. (2001) indicate a true GRB event 
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rate of 1 per year within D = lOOMpc. GRB980425/SN1998bw (z = 0.0085) is consistent 
with this true event rate. For anisotropic but geometrically standard sources, GRB980425 
is hereby not anomalous, but consistent with the trend shown in Fig. 4. 

The authors greatfully acknowledge input from Kevin Hurley. MVP thanks R. Mochkovitch 
for discussions on GRB980425. This research is supported by the LIGO Observatories, con- 
structed by Caltech and MIT with funding from NSF under cooperative agreement PHY 
9210038. The LIGO Laboratory operates under cooperative agreement PHY-0107417. This 
paper has been assigned LIGO document number LIGO-P030022-00-D. 
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Table 1. A sample of 33 GRBs with individually determined redshifts a 



GRB 


Rcdshift z 


Photon flux b 


Luminosity c 




Instrument 


(170000 

y fuzz® 


0.695 


10 


2.13 X 


10 5S 




Q A V l\\TT?r^ 

bAA/ Wrt 


y /UoUo 


U.ooO 


u.y ( 


3.24 X 


10 57 




oAA/ \ v r L 


n7nooo 

y fuozo 


0.9578 


1.5 


7.04 x 


10 57 


0.072 


"D V'TIP /AQA/T 

Ka Irj/ AbM 


y i 1214 


3.42 


1.96 


2.08 x 


10 59 


> 0.056 


Q A V /WTT?r^ 

bAA/Wrt 


nori/l oe; 
yoU4zo 


n unfit; 


n oft 


1.54 x 


10 53 




oAA/ W r L- 


yoUDlo 


1 nnfi 
l.Uyb 


U.O 


3.28 x 


10 57 


> U.lz/ 


QAV iWTT?C* 
oAA/ \ V r L- 


yoU /Uo 


u.yoo 


2.4U 


1.15 x 


10 5S 


U. loO 


HyV irj/ Aolvl 


nnm oq 

yyuizo 


1.0 


1 P. A 1 

ID. 41 


2.74 x 


10 59 


U.UoU 


oAA/ Wrt 


yyuoUo 


1.3 


18.56 


1.85 x 


10 59 




"D AT 1 /Dr 1 A 

r>Al / rUA 


yyuoiu 


i.oiy 


o. ID 


1.40 x 


10 59 


U.UOO 


oAA/ Vvrb 


yyo705 


0.86 








0.054 


bAA/ Wl O 


yyu7i2 


0.434 


11.64 


7.y7 x 


10 57 


> 0.411 


bAA/Wr O 


n o 1 riAO 

yyiius 


0.706 


11.2* 


2.48 x 


10 5S 


< 0.079 


TTl l\Wf~\ /TVT"C 

Uly/ KO/JNlL, 


yyiiio 


1.02 


67.5 


3.70 x 


10 59 


0.051 


BA1 /PCA 


r\r\r\ 101 

000131 


4.5 


1.5* 


3.05 x 


10 59 


< 0.047 


TTl l~\SC\ /MTI 1 

Uly/KU/JNhj 


000210 


0.846 


29.9 


1.03 x 


10 59 




bAA/WlO 


000301C 


0.42 


1.32* 


8.37 x 


10 56 


0.105 


ASM/Uly 


000214 


2.03 










SAX/WFC 


000418 


1.118 


3.3* 


2.27 x 


10 58 


0.198 


Uly/KO/NE 


000911 


1.058 


2.86 


1.72 x 


10 5S 




Uly/KO/NE 


000926 


2.066 


10* 


3.13 x 


10 59 


0.051 


Uly/KO/NE 


010222 


1.477 










SAX/WFC 


010921 


0.45 










HE/Uly/SAX 


011121 


0.36 


15.04* 


6.63 x 


10 57 




SAX/WFC 


011211 


2.14 










SAX/WFC 


020405 


o.6y 


7.52* 


1.58 x 


10 58 




Uly/MO/SAX 


020813 


1.25 


y.02* 


8.19 x 


10 5S 




HETE 


021004 


2.3 










HETE 


021211 


1.01 










HETE 


030226 


1.98 


0.48* 


1.35 x 


10 5S 




HETE 


030323 


3.37 


0.0048* 


4.91 x 


10 56 




HETE 


030328 


1.52 


2.93* 


4.31 x 


10 58 




HETE 


030329 


0.168 


0.0009* 


7.03 x 


10 52 




HETE 



a Compiled from S. Barthclmy's IPN redshifts and fluxes 
(http://gcn.gsfc.nasa.gov/gcn/) and J.C. Greiner's catalogue on GRBs local- 
ized with WFC (BeppoSax), BATSE/RXTE or ASM/RXTE, IPN, HETE-II or 
INTEGRAL (http://www.mpe.mpg.de/ jcg/grbgeb.html) 

b in cm _2 s _1 

c Photon luminosities in s _1 derived from the measured redshifts and observed 
gamma-ray fluxes for the cosmological model described in §2 

d Opening angles 9j in the GRB-emissions refer to the sample listed in Table I of 
Frail ct al.(2001). 

'Extrapolated to the BATSE energy range 50 - 300 keV using the formula given 
in Appendix B of Sethi et Bhargavi (2001) 
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Figure Captions 



Figure 1. Shown are three redshift distributions: the observed sample derived from Table 
1 (white), the true sample assuming the GRB event rate is locked to the star- formation 
rate (hachured), and the sample of detectable GRBs predicted by our model according to a 
log- normal peak-luminosity distribution function (grey). The continuous line represents the 
cosmic star formation rate according to a A— dominated cold dark matter universe. 

Figure 2. Shown is a simulation of the redshift distribution of the observed GRBs as a 
function of flux-limit, corresponding to various instruments including the upcoming SWIFT 
mission. The results are derived assuming the GRB event rate to be locked to the star- 
formation rate, using the best fit log-normal peak-luminosity distribution function used in 
Fig. 2. HETE-II tresholds are 0.21 (SXC), 0.07 (WXM) and 0.3 (FREGATE) in units of 

— 2 — 1 

cm s . 

Figure 3. Shown is a comparison between the flux distributions derived from our model 
(in grey) and from a sample of 67 GRBs observed with IPN (in wite). The fit serves as a 
test for our model assumptions, namely a log-normal GRB-luminosity function and a GRB 
event rate locked to the star-formation rate. 

Figure 4. Shown is a plot of the opening angle 9j of GRB-emissions versus redshift z in 
the sample of Frail et al. (2001), as derived from achromatic breaks in the GRB light curves. 
These results indicate an anticorrelation between 9j and z. For standard GRB-energies, 
this introduces a peak-luminosity function of GRBs which is correlated with the beaming 
factor 1/fb- This allows the beaming factor to be determined also in terms of the ratio 
of the unseen-but-true GRB event rate to the observed GRB event rate, using the current 
flux-limited sample of 33 GRBs with individually measured redshifts. 



unseen true probability 



flim = 0.04 ph/cm/s (SWIFT) 
flim = 0.2 ph/cm 2 /s (BATSE) 
flim = 1 ph/cm 2 /s 
flim = 10 ph/cm 2 /s (ULYSSES) 
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